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1 Introduction 



About twenty years ago an interesting class of the so-called quasi-exactly solvable (QES) 
potentials for which a finite number of eigenstates is analytically known was introduced 



Nowadays the QES problems attract much attention 0-19]. Several methods for 



generation of QES potentials have been worked out and as a result many QES potentials 
have been established. For example three different methods based respectively on a poly- 
nomial anzatz for wave functions, point canonical transformation and supersymmetric 
(SUSY) quantum mechanics are described in the paper |L2] . 



The SUSY is very useful tool for study of exactly solvable potentials. Note the papers 
|2(],|21| where SUSY procedure for constructing Hamiltonians either with identical spectra 



or with identical spectra, apart from a missing ground state, was given. This procedure 
may be repeated again and again to generate hierarchies of Hamiltonians whose spectra 
are related to each other. Some recent papers on this subject one can find in [f^-25|. For 
review of SUSY quantum mechanics see [2E,|27] 



At the first time the SUSY method for constructing QES potentials was used in ]|TC|-|T2|1 . 
The idea of this method is the following. Starting from some initial QES potential with 
n + 1 known eigenstates and using the properties of the unbroken SUSY one obtains the 
supersymmetric partner potential which is a new QES one with n known eigenstates. 

In our previous paper we proposed a new SUSY method for generating QES 
potentials with explicitly known two eigenstates. This method in contrast to the the 
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papers [p~0| p~2|] does not require the knowledge of the initial QES potential for generating 
of a new QES one. In the present paper we develop this SUSY method for constructing 
QES potentials with explicitly known three eigenstates. 



2 Witten's model of SUSY quantum mechanics 

Let us first take a look at the Witten's model of SUSY quantum mechanics. The algebra 
of SUSY in this case satisfies the following permutation relations 

{Q + 1 Q-} = H, (1) 
[Q±,H] = 0, 

(3^ = 0, 



where the supercharges read 



Q + = B-a + , Q- = B + a-, (2) 



are the Pauli matrices, 



W(x) is the superpotential. The Hamiltonian consists of a pair of standard Schrodinger 
operators H± 

where 

H ± = FB± = -±£ i + V ± (x), (5) 
V±(x) are the so-called SUSY partner potentials 



V ± (x) = - { W 2 (x) ± W'(xj) , W'ix) = ( 6 ) 
2 v ' dx 

Consider the equation for the energy spectrum 

H ± ^(x) = E±^(x), n = 0,l,2,.... (7) 

The Hamiltonians H + and H_ have the same energy spectrum except the zero energy 
ground state which exists in the case of the unbroken SUSY. This leads to two-fold 
degeneracy of the energy spectrum of H except for the unique zero energy ground state. 
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Only one of the Hamiltonians H± has the zero energy eigenvalue. We shall use the 
convention that the zero energy eigenstate belongs to H_. Due to the factorization of the 
Hamiltonians H± (see (§)) the ground state for H_ satisfies the equation 



B-^(x) = (8) 

with the solution 

^~[ x ) = C exp (- J W(x)dx\ , (9) 
here C is the normalization constant. From the normalization condition it follows that 

Bign(W(z)) = ±1, (10) 

when x — *> ±oo. 

The eigenvalues and eigenfunctions of the Hamiltonians H+ and H- are related by the 
SUSY transformations 

E~ +1 = E+, £ ~ = 0, (11) 

^- +1 (x) = -Ll^+(x), (12) 

' Ej, 



1>+{x) = fl-^Or). (13) 



J n+1 



The two properties of the unbroken SUSY quantum mechanics, namely, a two-fold 
degeneracy of the spectrum and the existence of the zero energy ground state are used for 
the exact calculation of the energy spectrum and wave functions (see reviews p6|, f2~7]j ) . 



3 SUSY constructing QES potentials 

We shall study the Hamiltonian H_ the ground state of which is given by (Q). Let us 
consider the SUSY partner of H_, i.e. the Hamiltonian H + . If we calculate the ground 
state of H + we immediately find the first excited state of H- using the degeneracy of 
the spectrum of the SUSY Hamiltonian and transformations (0), flT2]) , flilf). In order to 
calculate the ground state of H + let us rewrite it in the following form 

H + = H m + e, e > 0, (14) 

where 

H^=B+Bi, (15) 
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e is the energy of the ground state of H + since H W has zero energy ground state. 

As we see from (0), (0) the ground state wave function of H + is also the ground 
state wave function of and it satisfies the equation 

B^+(x) = 0. (16) 

The solution of this equation is 

$+{x) = C exp (- J Wi(x)dx^) . (17) 



Then using SUSY transformation fllT^D we can easily calculate the wave function of the 
first excited state of H_. Repeating the described procedure for H_ we obtain the second 
excited state for H . Continuing this procedure N times we obtain iV excited states. This 



procedure is well known in SUSY quantum mechanics p3| ,f21"|1 (see also reviews |26|. p7|) 



The wave functions and corresponding energy levels read 

^~(x) = C-B+...B+_ 2 B+_ 1 e X p^-Jw n {x)dx^ (18) 

n-l 

K = j>, ( 19 ) 

i=0 

where n = 1, 2, . . . , N. In our notations e = e, Bq = B ± , W (x) = W(x). Operators B^ 
are given by (|^) with the superpotentials W n (x). The equation (|TJ) rewritten for steps 

R{ n) = H (n+i) n — 0,1, . . . , N — 1 (20) 

leads to the set of equations for superpotentials 

W 2 n {x) + W&x) = ~ W^x) + 2e n , n = 0,1,..., N-l. (21) 

Previously this set of equations for W n (x) was solved in the special cases of the so-called 
shape-invariant potentials and self-similar potentials for arbitrary iV (see review plf). 
For N = 1 one can easily obtain a general solution of ( pip without restricting ourselves 
to shape-invariant or self-similar potentials. This solution was obtained in [[H]] in the 
context of parasupersymmetric quantum mechanics. 

In our recent paper |28| we construct nonsingular solution of (|2~lD for iV = 1 in order 



to obtain nonsingular QES potentials with two known eigenstates. 

In present paper we use the method proposed in |28| to solve the set of equations (|2l|) 



for iV = 2. It give us the possibility to obtain the general expression for wave functions 
and energy levels of QES potentials with three explicitly known eigenstates. 
Let us introduce new functions 

Wi n \x) = W n+1 (x) + W n (x), (22) 
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Wr\x) = W n+1 (x)-W n (x), n = 0,l,...,N-l. 

Then equations ([21]) read 

W' + {n \x) = W { _ n \x)W { +\x) + 2e n . (23) 

One can easily solve these equations with respect to the functions WA{x) obtaining the 
following expressions for superpotentials 

W n ( X ) = i(w^\ X )-^y-^\, (24) 



1 („ <n)t , , W^'(x)-2e n 



W n+1 (x) = - Wi n, (x) + 



2 I Wi n) (x) 



which lead obviously to the following set of equations for functions W+ 1 (x 



wrHx) + = ^ + (25) 



n = 0,l,...,iV-2. 



Thus the set of N equations (pi]) is reduced to the set of N — 1 equations fl25|). In 
the simplest case N — 1 equation (|25|) is absent and relations (|24D express just a general 
solution of equation (pip . To obtain a general solution of set ( pT|) in the case iV = 2 we 
have to solve one equation: 

rir , , W'(x)-2e ~ , . WL(x)-2ei . . 

where we have introduced for short the notation 

W + (x) = Wi°\x), W + (x) = w£\x). (27) 
It is easily to rewrite this equation as follows 

W+(x)W+(x)(W+(x) - W+(x)) - (W+(x)W + (x))' + 2(eiW+(x) + eW+(x)) = 0, (28) 



or 



u ( x ) I 7rr~T~\ - w +( x )) ~ u '( x ) + 2 I £xW+(x) + e V ^ x 



WAx) 




WAx 



(29) 



5 



where we have introduced the function 



U(x) = W+{x)W+{x). (30) 

We again came to the Riccati equation with respect to U(x). On the other hand it is 
an algebraic equation with respect to W + (x). Thus we can start from arbitrary function 
U(x) to construct the functions W + (x) and W + (x) which take the form 

_ 2U{x){U{x)+2e) 
W+[X} ~ U'(z)(l+K(x))' [31) 



W+(x) 

where 



U'(x) (1 + K(x)) 
2{U{x)+2e) 



TZ{x) = ±R(x), (32) 

R(X) : 



A U(x)(U(x) + 2e)(U(x)-2e 1 ) 



We mean square root R(x) as positively defined value, while function lZ(x) can be chosen 
in the form of R{x) or —R(x) within different intervals separated by zeros of function 
R(x). Note that just the possibility to choose different signs allows us, as will be shown 
in the section [5], to construct in a simple way new exactly solvable potentials using the 
known ones. 

Now we can obtain three consequent superpotentials W(x), W\{x) and W 2 (x) using 
the relations (pl|). Then using fllTf ) and ( |T9"D we obtain the energy levels and the wave 
functions of the first and the second excited states for if_ 



#r = e, Ez=e + e 1 , (34) 
^-( x ) = dW+(x) exp (- J W^dxj , (35) 

^-( z ) = C 2 ((W(a:) + W 2 (x))W+(x) - W' + ) exp (- j W 2 (x)dx 



The superpotentials W\{x) and ^(x) must satisfy the same condition (|T0|) as W(x). 
It leads to the same limitations for the functions W + (x) and W + (x). Both must be positive 
at infinity, negative at minus infinity and therefore each of them must posses at least one 
zero. Let us consider at first continuous superpotentials. As is seen from the fl24]) to avoid 
singularity of superpotentials the functions W + (x) and W + (x) ought have each only one 
zero |p8| 

W + (x ) = 0, W + (x ) = 
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at which they must satisfy the condition 



W^(x ) = 2e, W\.(x )=2e 1 . (36) 

Thus we have a number of limitations in the choice of function U(x) as a product of 
W + (x) and W + (x). There are two different possibilities of the choice of the function U(x). 
Either xq = xq that provides U(x) to have only one second-order zero point and to be 
positive at all the rest number line 

f U(x ) = 0, U'(x ) = 0, U"(x ) > 0, 

\ U(x) > 0, X ^ Xq, K ' 

or xq xq and therefore U (x) has two zero points and changes its sign as follows 
U(x) < 0, x G (min[xo, xq], max[xo, x~o\), 

U(xq) = U(x~o) = 0, U'(mm[xo,xo\) < 0, U' (max[xo, xq]) > 0, (38) 
U(x) > 0, x ^ [min[x , x ], max[x , x ]]- 

The sign of the function 1Z(x) in the expressions (|3l|) for the functions W+(x) and 
W+(x) should be chosen in a such way to ensure smoothness of these functions and an 
existence of one zero for each of them. A full analysis of the conditions which must satisfy 
the function U(x) to provide continuous superpotentials is rather boring and includes 
consideration of the behaviour of the functions W+(x) and W + (x) nearby possible zeroes 
of the expressions 

U(x)+2e, U(x)-2e u U'(x), R{x), (39) 

which is crucial for the continuity of the final superpotentials W(x), Wi(x) and Wi{x). 

We shall consider more closely the simplest case of functions U(x) which has one zero 
point and satisfies the conditions (pTf). The other condition which the function U(x) must 



satisfy is a consequence of the conditions (36) that connect the derivatives of the functions 



W + (x) and W + (x) with the energies e and 6\. One can easily find that 

U"(x ) = WI{xq)W + {xq) + 2W' + {xq)W' + {xq) + W + {xq)WI{xq) (40) 
= 2W^x )W[(x Q )=8ee l . 

The other obvious condition imposed on the function U (x) is positivity of the expression 
under the square root of the function R(x) 



U'{xf + 417(30(^(30 + 2e)(U(x) - 2eQ 

U'(x) 2 ~ { ' 

One can easily check that R(xq) = 0. Let us consider the case when point xq is unique 
zero of the function R(x). Then the only way to construct non-singular potentials is to 
chose function TZ(x) = R(x) over all the line in (|32|) . Besides we shall require the function 
R(x) to be smooth in the vicinity of point xq to avoid cusps of functions W + (x) and 
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W+(x) at this point. Note that such cusps would result in 5 like singularities of the final 
potential V^(x). 

Thus we obtain one of possible sets of conditions for the function U(x) allowing to 
construct non-singular QES potentials with three known eigenstates 

U(x) > V x ^ x , (42) 
U(x ) = 0, U'(x ) = 0, U"(x ) = 8eei, U'"(x ) = 0, 

U^(x ) = 64e ei ( ei - e), U^(x ) = 0, > 32(2e 2 - 13e Cl + 2e 2 ), 
R(x) > V x ^ x . 

Note that the most of exactly solvable potentials which are continuous over all the 
line satisfy these conditions. 



4 Examples 

One can easily check that the simplest functions U(x) yield us the well-known potentials. 
For example starting from U(x) = Aee\x 2 at e± = e we get harmonic oscillator potential. 
Another simple function U(x) = 4eei tanh 2 x at t\ = e — 1 leads us to the well-known 
exactly solvable Rosen-Morse potential. One more simple example U (x) = 4eei sinh 2 x at 
ti = e + | reproduces the special case of the well-known quasi exactly solvable Razavy 
potential [|J. 

Let us consider more complicated examples leading to new QES potentials. We shall 
start from the function 



U(x) = Aee lX 2 \±4^, (43) 



b 2 x 



where a and b are real parameters. 

Due to the conditions for U^(x ) and U^^Xq) from (021) we get 



and 



respectively; here 



a 2 = b 2 + 2 -{e 1 -e) (44) 



13 e 15 3 a— ei 13 e 15 3 i— . . 

VA< — < h-VA (45) 

4 b 2 8 8 ~ b 2 ~ 4 b 2 8 8 y ' 



e _/e^ 2 



A = 25 - 60— + 68 f - ) . (46) 
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Because a and b are real numbers it follows from ([44]) that 



---<-• (47) 

b 2 2~b 2 1 J 

It is easily to check that 

13 e 15 3 rr e 3 . 

VA < (48) 

4 b 2 8 8 - b 2 2 y 1 

and therefore (|45|) together with (f44|) lead to the following inequality 

e 3 ei 13 e 15 3 r— . , 



Thus in order to obtain nonsingular potentials the parameters of function ([43]) must 
satisfy conditions (fP|) and ([13). Besides we obviously must require the parameters e and 
ex to be positive. 

We shall omit the general expression for the potential VL(a?) as it is huge and rather 
useless. It is easily to show that there are only two sets of positive e and e\ allowing to 
resolve the root in the function R(x) (Cases 1 and 2) and therefore to simplify significantly 
the final expressions. The other simplified expression (Case 3) we shall obtain putting 
a = that corresponds to the lowest value of the parameter e% in inequality ([49|). 

4.1 Case 1 

In the case 

e = h\ ei =(l + Vs)b 2 (50) 



2 ' V2 

square root in R(x) can be resolved and we obtain the potential 

V-(x) = (-?-2V3 + 3(I + ^f.V+ (51) 

7V3-12 w 



2-V3 + b 2 x 2 (2 - + 62 x 2 

and the eigenfunctions 

M*) = Ce-^^\2-^ + b^)^ 



(52) 



(3+2v^), 2 2 / _ oX i , ^- 1 ) 

^(x) = c ie -^-^ bx (2-V3 + b 2 x 2 ) 2 x, (53) 
V*(a:) = C 2 e"^^ 2 (2 - v/3 + 6 2 x 2 ) ^ f 2 - V3 - &V) (54) 
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corresponding to the three lowest levels 



Eq=0, E^=h\ E 2 =(3 + V3)b 2 . (55) 



2 

To simplify these expressions let us make the substitution 

b 2 = (2 - Vz) c 2 . (56) 

Thus we get 



, /3 2 2 3- V3 2^-3 7^/3 3\ 

V_m = c \ -c x H r-^r H 9 h- , (57) 

U l + c 2 x 2 l + c 2 x 2 ) 2 4 2 ' 



£- = 0, £f = 3 f 1 - ^ c 2 , £? 2 " = (3 - V3) c 2 , 



(5f 



^-(x) = ^-^(l+cV)^, (59) 



$-(x) = Cje * cV (l + A 2 ) 2 (60) 
^-(x) = C 2 e-^ 2 (l + cV) (l-cV). (61) 



It is worth to note that the obtained potential fl5"T|) is double-well one. As well known, 
double-well potentials has been used extensively to model a wide range of natural phe- 
nomena. 

4.2 Case 2 

The other set of positive e and 6\ resolving the root in the function R(x) reads 

e = h 2 , £1 = h 2 . (62) 



It leads to the well-known supersymmetric partner of a harmonic oscillator [22,23.33 



1 . 36 2 &V 46 2 26 2 . , 

V-[x) = — + — ;9 9 , 9 + - ■ , 9 9 . (63) 

4 8 (1 + 6 2 x 2 ) 2 1 + 6 2 x 2 

This potential is exactly solvable although we have got it using the procedure for con- 
structing QES potentials with three exactly known eigenstates. We will also obtain the 
same potential in the section ||] using the other procedure allowing to construct new exactly 
solvable potentials (see flSiP). 
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4.3 Case 3 



Let us consider another particular case of the potential under consideration which is rather 
simple. We put a = to reduce the function U(x) ([43]) to the form 



U(x) = 4eei— (64) 
1 + trar 

Due to the relation fl44]) we obtain the connection between the energies e and E\ 

e = e 1 + \b\ (65) 



Note that in accordance with relation (f49|) such a choice of parameters corresponds to the 
minimal value of t\ at given e and b. It leads directly to the potential 

V-{x) (a 2 + 3a + l) 2 1 2 2 

& 2 8a(a + l) + l + 6 2 x 2 (1 + b 2 x 2 ) 2 p(x)(l + b 2 x 2 ) 2 1 J 

a 2 + a-l (a 2 + ct + l) 3 (a 2 + a + l) 2 



+ 6 2 x 2 ) 8a(a + l)p 2 (x) Ap 3 (x) 
and the following eigenfunctions 

*T(*) = Cbg±fe-*^( 1+ ^), (67) 

^-(a:) = d— ^ e -3K*)(i-|+5^) (68) 

- 1 

^-(ar) = C 2 p{x) ~° l ~ l e-^-^M, (69) 
p(x) - 1 

corresponding to the levels 

E- = 0, E ^ = e 1 + h 2 , E 2 =2 ei + h 2 , (70) 

where 



p(x) = v /l + a(a + l)(l + 6 2 x 2 ), a = 1 + 2^. 



Obtained potential (p6|) has one minimum at x = and tends to constant for x — > ±oo. 
For the case a > 1 all the obtained wave functions (|67|-|6"9"|) are square integrable and there 
exist at least three bound state in the well. For (V5 — l)/2 < a < 1 only two lower wave 
functions (p7|),(pB|) are square integrable and we have just two bound states, while for 
ol < (\/5 — l)/2 only the ground state wave function d67|) remains square integrable and 
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the potential has only one bound state. 



5 Constructing exactly solvable potentials 

Although we have developed our scheme to construct new QES potentials it seems to be 
also of use for constructing SUSY partner potentials for the exactly solvable ones. Let us 
start from some exactly solvable potential V-(x) 



Vjx) 



i (\V\x) - W'{x)) , (71) 



for which we know three first superpotentials W(x), W\(x) and W^(^) satisfying set 
of equations fl2T|). One can easily construct the functions W + (x) = W(x) + Wi(x), 
W+(x) = W\(x) + W 2 (x) and U(x) = W+(x)W + (x). Substituting function U(x) into the 
expressions (|31| ) we can obviously reproduce the functions W + (x) and W+(x) choosing 
corresponding signs in the function ([32]) which we shall denote as TZq(x). Beside the 
functions W + (x) and W+(x) we can obviously obtain another pair of functions W + (x) 
and W+(x) given by the same expressions ([31]) with the only difference that here we 
choose the function (|3"2"D with the opposite sign to TZq(x), 



TZ(x) = -Ko(x) (72) 
The new functions W + (x) and W+(x) satisfy the same equation (2H) and just they allow 



us to construct new exactly solvable potentials. 

Let us consider these functions in more detail. It is easily to show that both functions 
W+(x) and \V + (x) are negative at the infinity and positive at the minus infinity. Ex- 
plicit calculations show that it provides the same behaviour of the superpotentials W(x), 
Wi(x) and W 2 {x) which can be obtained by substitution functions W + (x) and W + (x) 
into relations (j2"4|). We shall omit explicit expressions for the superpotentials Wi(x) and 
W2(x) noting only that both of them are singular while the superpotential 

_ U'(x)(n (x)-1) 8ee 1 -U"(x) + 2U(x)(A(e 1 -e)-3U(x)) 
VV[X) ~ 2(2e + U(x)) + 2U'(x)n (x) { ^ } 

has no singularities if only the initial potential V_(x) is non-singular. Using superpo- 
tential W(x) we can construct in a standard way the pair of the Hamiltonians and 
7i + . Their properties will be very similar to that of the Hamiltonians H_ and H + with 
the only difference that now Hamiltonian 7i + will have zero-energy ground state with 
corresponding eigenfunction 

ip+(x) = Cexp (J W(x)dx ) j . (74) 
All the higher eigenvalues of Hamiltonians 7i + and 7i- will coincide and corresponding 
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eigenf unctions will be connected as follows 



"n+1 

Let us consider few the simplest explicit examples. In the case of a harmonic oscillator 
for which all the superpotentials and distances between the energy levels read 

W n (x) = ex, e n = e (76) 

we get corresponding function 

U(x) = 4eV (77) 

and we obtain 

Ho{x) = R(x) = 2ex 2 . (78) 
Choosing now 1Z(x) = —TZq(x) and using (fT3|) we get 

K _|_ Of T 2 

W(x) = -ex--—;, (79) 



1 + 2ex 2 

leading to the following SUSY partner potentials 



V-(a?) = —+2", (80) 
, . e 2 x 2 4e 8e 3e , > 

v + w = — + IT2^ - ITS?? + 2' (81) 



Thus we have obtained the exactly solvable potential ( |STD which is a SUSY partner to the 
harmonic oscillator potential (BO). This potential (|8T| ) after substituting e = 6 2 /2 coincides 
with potential (|63"D obtained in the section |[ Let us remind that the potential V+(x) has 
zero-energy level with corresponding wave function (|74]) . Therefore we can now treat 
the potentials V_(x) and V+(x) as U + (x) and V-(x) respectively, that is V + (x) = V_(x) 
and U_(x) = V+(x). The superpotential corresponding to the pair V-(x) and V + (x) is 
= — W(x). Because the upper SUSY partner V+(x) is a harmonic oscillator we can 
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easily build up all the hierarchy of superpotentials satisfying equations (^Tj) : 

5 _i_ 2ex 2 

W (x) = ex———, W n (x) = ex, n = l,2,... (82) 
1 + lex 1 

The corresponding distances between the energy levels read 

e = 3e, e n = e, n = l,2,... (83) 

The obtained potential (|8lD can be used for constructing another exactly solvable 
potential. Starting from the first three superpotentials fl8"2"|) we obtain in the same way as 
before the next pair of potentials 

V-{x) = —+2", (84) 

2^2 o^r. 2 q\ QS/1^2^,2 



ex 



8e(2ex 2 -3) 384eV 7e 



2 + 3 + 12ex 2 + 4e 2 x 4 + (3 + 12ex 2 + 4e 2 x 4 ) 2 + 2 ' ^ 

Repeating this procedure many times we obtain the following pairs of SUSY partner 
potentials 

V + (n,x)=V4n,x) = ^ + (in - ~) e. (86) 

V4n,x)=V + (n,x) = ?f + 8en(2n - 1) H ^~ff^ (87) 

2 H 2n {i^/ex) 

Wcn 2 ( H2n - l{iV ~ ex) X I (4n ~ 1)£ 



H 2n {i\fex) J 2 



where H n (x) - - Hermite polynomial. Note that the case n = 1 in the potential fl87Q 
corresponds to (^1|) and the case n = 2 corresponds to (|S5D. The potentials V-(n, x) (p7| ) 
are just the special cases of the SUSY partner potential of a harmonic oscillator obtained 
by Sukumar |^T[ and they were previously obtained by Bagrov and Samsonov |22|, f23| via 



the Darboux method and latter by Junker and Roy |33| within the SUSY approach. 



The application of the same procedure for the Morse and Rosen-Morse potentials 
as well as for radial harmonic oscillator and hydrogen atom provides chains of exactly 
solvable potentials. All these potentials are just special cases of the potentials obtained 
in by Junker and Roy. Nevertheless they seem to be interesting because all of them 



are expressed in terms of the elementary functions only. For example starting from the 
Morse potential for which corresponding superpotentials and distances between the energy 
levels read 

W n (x) = e + 1/2 - n - e~ x , e n = e - n, (88) 
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we construct the function 



U(x) =4(e-e- x )(e-l-e~ x ). (89) 
An explicit calculations show us that we should take function TZq(x) in the form 

^ . . 4e~ 2x + 6(1 - 2e)e~ x + 3(1 + 4e(e - 1)) - 2e(l - 3e + 2e 2 )e x , nn , 

■io(X) = ; r (9U) 

v ' 2e~ x + (1 - 2e) ( J 

to reproduce Morse potential. Function TZq(x) does not coincide with R{x) in this case, 
the former is negative within the interval limited by zeros of the expression ( |90|) . Following 
the same procedure as in the case of harmonic oscillator we obtain such a sequence of 
potentials 

v _ {0 , x) = ii±^ + ^- 2 '; +1 ^ , (si, 

y_ (l , x) = <l±if>! + e - - % - *~ ,92) 



+ 



2(e(2e - l)e x - 2) 8((2e - l)e x - 1] 



e(2 - 2(2e - l)e x + e(2e - l)e 2:r ) e(2 - 2(2e - l)e x + e(2e - l)e 



2xA2 ■ 



= (l+*g + - - 3)e- (93) 



+ 



+ 



(2e - 3) (8e x - (e - 1) (48e 2a; - (2e - 1) (36e 3a; - 16ee 4x ))) 
4 - (2e - 3) (8e* - (e - 1) (12e 2 * - (2e - 1) (4e 3 * - ee 4 *))) 

((2e - 3) (8e x - (e - 1) (24e 2x - (2e - 1) (12e 3x - eAe ix )))f 
(4 - (2e - 3) (8e* - (e - 1) (I2e 2x - (2e - 1) (Ae 3x - ee 4 *)))) 2 ' 



Corresponding SUSY partners are following Morse potentials 

xr / \ (l + 2e) 2 e~ 2x - 2(e - 2n)e~ x , nA . 
V+(n,x) = g + ^ ^ . (94) 

We can proceed with this procedure as long as necessary. 

The most interesting fact is that at each step of the suggested procedure for the all 
mentioned above potentials 

H- = #[ 2) + e + ei . (95) 
One can easily check that due to the connection ( O ) and (p0|) 

E n )+ = E nll = En+1 - e l = E n+2 ~ e l _ 6 = E n+3 ~ 6 l ~ C " ( 96 ) 

and therefore the energy levels of the Hamiltonian 7i_ coincide with that of the Hamilto- 



15 



nian H- saving the three lowest levels of the latter which are not present in the spectrum 
of the former 



£- = E~ +3 , n=0,l,... (97) 

It gives us immediately the energy spectrum of the new Hamiltonian Ti, + 

£ =0, S- = E~ +2J n = l,2,... (98) 

Besides the connection (|95| ) allows us to obtain easily all the eigenfunctions of the excited 
states for the new exactly solvable potential TC + . Using relation (|75|) we obtain 

cp+(x) = C+ (Jt + W(x)j ^lt(x), n = 1, 2, . . . (99) 

Thus for all mentioned exactly solvable potentials we can construct a sequence of 
Hamiltonians H n = —\^-% + V(n,x) with the same energy levels as H except the 2n 
lowest excited states of the latter. All of them possess zero-energy ground state. 



6 Conclusions 

We have obtained a general solution of the set of equations ([H]) for N = 2 given by 
the expressions ( pi]) and (Q) (n = 0,1). Thus we can write down explicit expres- 
sions for superpotentials Wq(x), W\{x) and W%{x) and afterwards for potential vl°\x) = 
(Wq(x) — Wq(x))/2 being just a general expression for a QES potential with three explic- 
itly known eigenstates. General expressions for the corresponding eigenfunctions are also 
presented (||), ([35|). The QES potential is expressed in terms of the distances e, ei between 
neighbouring energy levels and an arbitrary function U(x). To ensure nonsingularity of 
the potential we need to put a number of limitations on the function U(x). Using this 
expression we have obtained some new QES potentials. In special cases our potentials 
reproduce those studied earlier. 

There obviously arises a question if the suggested scheme could be generalized to 
construct QES potentials with the number of explicitly known eigenstates large than 
three. In such a case we have the set of equations fl25|) consisting of N — 2 equations. 
Let us remind that N is number of excited states of QES potential which we would like 
to construct. In order to reduce this set of equations we can proceed with the scheme 
described in section |3] (eqs. (PT|)-(PT|)).The solution of each equation of the set (|23| ) can be 



written down in the form (R|) , where W + (x) is replaced by VU| n) (x), W+ (x) is replaced by 
Wi n+1 \x) and U n (x) = wf{x)W [ + n +1) (x) is instead of U(x). Two neighbouring equations 
of the set (|25|) yield us two different expressions for the same function in terms of the 
functions U n (x) and U n+ i(x) correspondingly. Thus we obtain following set of equations 
for functions U n (x) 



U' n {x) (1 + K n {x)) 2U n+1 {x){U n+1 {x) + 2e n+1 ) 



2(U n (x) + 2e n ) U' n+1 {x) (1 + K n+l (x)) 
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n = 0,...,JV-3, (100) 



where 

^M.xL A U n (x)(U n (x)+2e n )(U n (x)-2e n+1 ) 

MI] ( + U' n {xf ' 1 J 

Number of equations in this set is one less than number of equation in the set (|25|) 
for functions W+ 1 (x) and correspondingly it is two less than the number of the initial 



equations fl2l|) for surperpotentials. In the case N = 3 we have just one equation which 
we need to solve to construct QES potentials with four known eigenstates. Thus one 
can see that the obtaining of general expression for QES potentials with more than three 
explicitly known eigenstates is essentially more complicated. 

Another point is that the suggested scheme allows one to construct in a simple way the 
sequences of SUSY partner potentials of exactly solvable ones. At each step we obtain a 
new exactly solvable potential with identical spectrum, apart from missing two the lowest 



excited states. Note that this approach in contrast to method used in p0|-p3|, |32| , |33]j 
does not require knowledge of general solution of corresponding Schrodinger equation for 
initial potential. It would be interesting to apply the same approach to the known QES 
potentials with explicitly known N eigenstates to construct new QES potentials with 
picked out two the lowest excited states. But the latter case is more complicated than 
the case of shape-invariant potentials and it will be the subject of a separate paper. 
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